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Abstract. I review a proposed pattern of the light scalar mesons with qq mesons and glueball above 
1 GeV and tetraquark mesoniums below 1 GeV. Several challenges and caveats of calculating these 
light scalar mesons with dynamical fermions are discussed. 
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"V INTRODUCTION 

The pseudoscalar, vector, axial, and tensor mesons with light quarks (i.e. w, d and s) 
are reasonably well known in terms of their SU (3) classification and quark content. 
C\) ■ The scalar meson sector, on the other hand, is much less understood in this regard. 

There are 19 experimental states below 1.8 GeV which are more than twice the usual 
qq nonet in other sectors. We show in Fig. Q] the experimentally known scalars includ- 
ing c(600), ?c( 800), and /o(1710) which are better established experimentally nowa- 
days 0], 0]. The recent theoretical advance [d] in identifying a(600) as a %% resonance 
by solving the Roy equation has settled the question about the existence of c(600). 
Nevertheless, there are still a number of puzzling features regarding the ordering of 
ao(1450) and £q(1430) with respect to their counterparts in the axial-vector and tensor 
sectors, the narrowness of ao(980) and /o(980) in contrast to the broadness of a(600) 
and k(800), etc |4Q. We shall first review a emerging pattern of the scalar mesons below 
1.8 GeV based on quenched lattice calculation and phenomenology and then discuss the 
challenges and caveats of full QCD calculation of these scalar mesons on the lattice. 



PATTERN OF LIGHT SCALAR MESONS 

The unsettling features regarding the nature of an(1450) and ^(1430) are tentatively 
resolved in a recent quenched lattice calculation pD with overlap fermions for a range 
of pion masses with the lowest one at 180 MeV When the quenched ghost states, which 
correspond to nr\ and nr\' scattering states in the dynamical fermion case are removed, 
it is found that ciq is fairly independent of the quark mass. In other words, below the 
strange quark mass, ao is very flat and approaches ao(1450) in the chiral limit. This 
suggests that SU (3) is a much better symmetry in the scalar meson sector than the 
other meson sectors and that both ao(1450) and ^(1430) are qq states. Furthermore, 
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FIGURE 1. Spectrum of scalar mesons together with it, p , a\ and a%. 



/o(1500), by virtue of the fact that it is close by, should be a fairly pure ST/ (3) octet 
state, i.e. f octet = { u u + dd — 2ss) / y/6. 

Based on the lattice findings, a mixing scheme for the isoscalar /o( 1370), /o( 1500) 
and /o(1710) - a glueball candidate, with slight SU(3) breaking was developed and 
successfully fit to the decays of pseudoscalar meson pairs as well as various decays 
from J/^ [6|]. Some of the robust and conspicuous features of this mixing scheme are 
the following: 

• /o(1500) is indeed a fairly pure octet with very little mixing with the flavor singlet 
and the glueball. /o(1710) and /o(1370) are dominated by the glueball and the qq 
singlet respectively, with ~ 10% mixing between the two. This is consistent with 
the experimental result r(//Y -»• 7/0 (1710)) ~ 5T(J /m -> 7/0 (1500)) H which 
favors /o(1710) to have a larger glueball content. 

• The ratio r(/ (1500) -> ^Z)/r(/ (1500) -> %%) = 0.246 ±0.026 is one of the 
best experimentally determined decay ratios for these mesons [1]. If /o(1500) is a 
glueball (i.e. a flavor singlet) or ss, the ratio will be 0.84 or larger then unity. Either 
one is much larger than the experimental result. On the other hand, if /o(1500) 
is foctet, then the ratio is 0.21 which is very close to the experimental value. This 
further demonstrates that /o(1500) is mainly an octet and its experimental decay 
ratio can be well described with a small SU (3) breaking 

• Because the nn content is more copious than the ss in /o(1710) in this mixing 
scheme, the prediction of T(J/x\r -> Q)/o(1710))/r(J/y 0/o(171O)) = 4.1 is 
naturally large and consistent with the observed value of 6.6 ± 2.7. This ratio is not 
easy to accommodate in a picture where the /o( 17 10) is dominated by ss. One may 
have to rely on a doubly OZI suppressed process to dominate over the singly OZI 
suppressed process to explain it D_ZD • 
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FIGURE 2. Pattern of light scalar mesons - a tetraquark mesonium nonet below 1 GeV, an almost pure 
SU(3) qq nonet and a nearly pure glueball above 1 GeV. 

The mesons below 1 GeV were suggested to be tetraquark mesoniums Q from the MIT 
bag model [8] and potential model [|9j,LL0|] studies. A recent lattice calculation [5] with 
the overlap fermion on 12 3 x 28 and 16 3 x 28 quenched lattices with the two-quark-two- 
antiquark interpolation field my^mmy^V has confirmed the existence of such low-lying 
scalar tetraquark mesonium at ~ 550 MeV. This strongly suggests that it is the a (600). 

Combining the lattice calculations of ciq ( 1 450), Kq( 1430) and a (600) and the mixing 
study of /o (1370), /o( 1500) and /q(1710), a classification of the scalar mesons be- 
low 1.8 GeV was proposed fl]. Those below 1 GeV, i.e. fj(600),a (980),/ (980) 
and K"(800) form a nonet of tetraquak mesoniums; those above 1 GeV, i.e. 
flo( 1450), ££(1430) and / (1500) form a fairly sure SU(3) octet; and / (1370) 
and /o(1710) are good ST/ (3) singlet and glueball respectively, with ~ 10% mixture 
between the two. 

We should stress that this is not finalized. It should be scrutinized in future exper- 
iments, such as high statistics J D, B decays and pp annihilations. Furthermore, 
most of the the lattice results which led to the above proposed pattern were based on 
quenched calculations. There are loose ends that need to be tightened, come dynamical 
fermion calculations. In the following, we shall enumerate a number of challenges and 
the associated caveats in calculations with dynamical fermion configurations. 

CHALLENGES AND CAVEATS OF FUTURE LATTICE 
CALCULATIONS WITH DYNAMICAL FERMIONS 

In the quenched lattice calculation of ao with light quarks which correspond to m K < 500 
MeV, the quenched %r\ ghost states are lower than the ao(1450) and, thus, dominate the 
long time behavior in the scalar correlator with a non-unitary negative tail. This has to be 
removed fiQjltll] before the physical qq( 1450) is revealed. These ghost states turn into 



These are two-quark and two-antiquark mesons which have been referred to as four-quark states, meson 
moleculars, mesoniums, and tetraquark states. We shall call them tetraquark mesoniums so as to avoid 
implication on the nature of possible spatial and color clustering. 



physical two meson scattering states in a full QCD calculation with the same valence and 
sea quark masses H This causes some difficulty in the fitting of scalar meson correlators 
and has been mentioned by S. Prelovsek [| 1 3M in this workshop. In the following, we 
shall point out several caveats and challenges facing the scalar meson calculation with 
light dynamical fermions. 



a (1450) and K*{ 1430) 
There are several N± — 2 dynamical fermion calculations of ao with the interpo- 



lation field II12L ll4L U51 [lq. uJtt. Save for Ref. u2\\ which, upon removing the partially 
quenched ghost %r\2 state, found ao to be at 1.51(19) GeV, the others |jl5Lll6L ll7ll found 
the lowest states at the chiral limit to be ~ 1 GeV, suggesting that ao(980) is the qq state. 
As pointed out in Ref. [4], this is most likely an untenable interpretation. If ao(980) is 
indeed a qq state, or has a sizable coupling to the interpolation field, then replac- 
ing the ujd quark in the ao interpolation field with s will place the corresponding su at 
~ 1100 MeV. This is far (i.e. 300 MeV) away from each of the two experimental states 
Kq (1430) and fc(800) (see Fig. 1). The likely resolution, we think, is that the state found 
at ~ 1 GeV is the nr\2 scattering state. Since T72, the r\{r\') in the two-flavor case is 
predicted to be ~ y / 2/3m JJ / = 782 MeV in the large N c analysis with U(l) anomaly, the 
weakly interacting %r\2 will be near the state seen at ~ 1 GeV. In other words, this %r\2 
scattering state is the dynamical fermion realization of the corresponding ghost state 
in the quenched approximation. Parallel to the lesson learned in lattice calculations of 
pentaquark baryons tl8fl, one has to include the multi-hadron states in addition to the 
physical resonances when fitting the two-point correlators for the excited spectrum. In 
the case of ao in the realistic Nf = 2 + 1 case, one needs to include %r\ , nr\' , in addition 
to the physical ao(980), and «o(1450). This can be achieved with the sequential empir- 
ical Bayes method for curve-fitting Ii22ll or the variational approach. Furthermore, one 
needs to distinguish the two-particle scattering states from the one-particle resonances. 
One way to distinguish a two-particle scattering state from a one-particle state is to ex- 
amine the 3 -volume dependence of the fitted spectral weight |fl9l I20I Bl . Another way 



is to impose a 'hydrid boundary condition' on the quark propagators 112 ill . No attempt 
has been made to identify the scattering Tzr\(r)') states so far. This has to be carried out 
before one can reasonably reveal the quark content of ao(980) and ao(1450). 



/o(980),/o(1370),/ (1500) and/ (1710) 

In addition to the complication of two-meson scattering states (in this case 
TZTZ,KK,r\r\,r\r\'), one needs to calculate the correlators with disconnected inser- 
tions (D.I.) in addition to the connected insertions (C.I.) as in the ao case. This is 
to reflect the fact that these isoscalar mesons have annihilation channels. The usual 



Otherwise, it is considered to be a partially quenched calculation. 



approach of adopting the noise HH] to estimate the quark loops in the disconnected 
insertion makes the calculation much more expensive than the connected insertion one. 
One caveat with the noise estimator is that the signal falls exponentially in the meson 
correlator; whereas, the variance of the noise estimator approaches a constant at certain 
time separation Q23[|. If one were to fit the time window where the variance of the noise 
levels off, the shoulder effect of the correlator could result in an unphysically light 
effective mass. In view of this, the very light mass from the D.I. part of the correlator 
in the fo calculations [[Til I2H] should be subjected to the examination as to whether it is 
the nn scattering state or due to the shoulder effect. 



Glueball 

The continuum and large volume limits of the quenched calculation places the scalar 
glueball at 1710(50)(80) MeV [25]. This is very close to the viable experimental glueball 
candidate /o(1710). To verify this in the full QCD calculation is, however, non-trivial. 
Whatever interpolation one adopts, one has to disentangle the glueball from all the 
lower-lying /o states and the KK,KK,r\r\ and r\r\' two-meson states. 



Tetraquark Mesoniums 

If the nonet below 1 GeV in Fig. 1 are indeed dominated by the q 2 q 2 tetraquark 
mesonoums, one can access them through the ^75 ^^Vy^V operator or other four-quark 
operators with the same quantum number. In the case of ao(980) and /o(980), the 
two-meson threshold, i.e. KK is close by. One may need a good variational method 
in order to disentangle them. By virtue of the fact that ao(980) and /o(980) are nearly 
degenerate, the D.I. should be small compared to the C.I. It should be confirmed in full 
QCD calculation. 

qq meson vs q 2 q 2 tetraquark mesonium 

The notion of qq or q 2 q 2 meson is primarily a quark model concept of the valence 
quark content. How does one distinguish them in lattice QCD with interpolation fields? 
So far, neither ao(980) nor a(600) is coupled to the x ¥ x ¥ interpolation field in the 
quenched approximation with discernable signal [5]. If this is not true in full QCD 
calculation with light dynamical fermions, this will complicate matters substantially. 
One will need both qq and q 2 q 2 types of operators with a large basis in the variational 
calculation to identify states and; moreover, to distinguish the one-particle states from 
the multi-meson scattering states. 



CONCLUSION 



We summarized the pattern of light scalar mesons emerged from quenched lattice cal- 
culations and the study of mixing and decays of /o (1370), /o( 1500) and /o(1710). We 
have discussed the subtlety and challenges of calculating them in full QCD with light 
dynamical fermions. In particular, if they couple strongly to both qq and q 2 q 2 types of 
interpolation operators, the interpretation of the underline pattern will be considerably 
more complex. We hope that Nature is only subtle but not malicious. 
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